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ABSTRACT 

Optical integral field spectroscopy (IFS) combined with Hubble Space Tele- 
scope (HST) WFPC imaging were used to characterize the central regions of 
the Seyfert 1 radio galaxy 3C 120. We carried out the analysis of the data, 
deriving intensity maps of different emission lines and the continua at different 
wavelengths from the observed spectra. Applying a 2D modeling to the HST 
images we decoupled the nucleus and the host galaxy, and analyzed the host 
morphology. The host is a highly distorted bulge dominated galaxy, rich in sub- 
structures. We developed a new technique to model the IFS data extending the 
2D modeling (3D modeling hereafter). Using this technique we separated the 
Seyfert nucleus and the host galaxy spectra, and derived a residual data cube 
with spectral and spatial information of the different structures in 3C 120. Three 
continuum-dominated structures (named A, B, and C) and other three extended 
emission line regions (EELRs, named Ei, E2 and E3) are found in 3C 120 which 
does not follow the general behavior of a bulge dominated galaxy. We also found 
shells in the central kpc that may be remnants of a past merging event in this 
galaxy. The origin of Ei is most probably due to the interaction of the radio-jet 
of 3C 120 with the intergalactic medium (Axon et al. 1989; Sanchez et al. 2004a). 
Structures A, B, and the shell at the southeast of the nucleus seem to correspond 
to a larger morphological clumpy structure that may be a tidal tail, consequence 
of the past merging event. We found a bright EELR (E2) in the innermost part 
of this tidal tail, nearby the nucleus, which shows a high ionization level. The 
kinematics of the E2 region and its connection to the tidal tail suggest that the 
tail has channeled gas from the outer regions to the center. 
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1. 



Introduction 



Radio galaxies are excellent laboratories to disentangle the role of central active nuclei in 
host galaxies evolution and the relation of intergalactic environment and the central activity. 
The complex structure of radio galaxies put severe troubles to carry out specific optical 
studies in this matter. Fortunately, the recent technical and computational improvements 
on integral field spectroscopy (IFS hereafter) open new frontiers to multitude of astronomical 
projects. In the particular case of radio galaxies, IFS data allow to separate the different 
components and to perform an independent analysis of the various subsystems and their 
evolution. 

3C 120 is a Seyfert 1 galaxy at a redshift of 2; = 0.033. It has been studied extensively 
at many wide-range wavelengths, with special attention in X-ray and radio frequencies (e.g. 
Walker et al. 1987; Moles et al. 1988; Walker 1997; Ogle et al. 2004). Ahhough 3C 120 
was classified as an early-type galaxy based on a visual inspection of ground-based images 
(Zwicky & Zwicky 1971), Sargent (1967) reported a faint spiral structure. However, its 
optical morphology is still not clear showing an elliptical shape with some peculiarities 
(Hansen, Jorgensen, & Norgaard- Nielsen 1987) and could be a merger remnant galaxy 
(Moles et al. 1988). Spectroscopic analysis has shown that 3C 120 presents a rotation curve 
most likely indicating the presence of an undetected disk (Moles et al. 1988). Based on this 
result, 3C 120 is normally quoted as an early-type spiral or SO galaxy (Moles et al. 1988). 

Although Seyfert galaxies are generally considered as radio-quiet objects, 3C 120 is 
an active radio-loud source showing a superluminal one-sided jet that extends ~25" out of 
the core (Readhead et al. 1979; Walker et al. 1988; Walker 1997). Different authors (e.g. 
Baldwin et al. 1980; Pcrez-Fournon et al. 1986; Soubeyran et al. 1989; Hua 1988) reported the 
existence of several continuum and emission line dominated structures in this object. Some 
of these structures have a clear relation with the radio-jet. Hjorth et al. (1995) detected 
a continuum dominated optical counterpart of the radio-jet using deep broad-band ground 
based images. On the other hand, an EELR (Ei) associated with the jet was found at ~5" 
west of the nucleus (Hua 1988; Axon et al. 1989; Sanchez et al. 2004a). The nature of the 
remaining structures is still unclear (e.g. Perez- Fournon et al. 1986; Moles et al. 1988). 

This paper compiles data from IFS optical observations of 3C 120 using the optical fibers 
system INTEGRAL covering a central region of ~ 16" x 12" field-of-view. We combined our 
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IFS data with high- resolution Hubble Space Telescope (HST) images. We describe the optical 
fibers instrument, observations, and data reduction in section §2. In section §3, we describe 
the emission line profiles in the circumnuclear region of 3C 120 and we present 2D spectra 
diagrams for different emission lines (appendix A). In the following sections, we present 
the stellar morphology and the structures detected (section §4.1), the average colors of the 
galaxy and the colors of those structures (section §4.2), and the emission- line morphology 
and ionization structure (section §4.3). In section §5, we apply a new technique to separate 
the spectra of different components coexisting in 3C 120 and analyze the properties of each 
component. We present the velocity field of the central region of 3C 120 and the discussion of 
several kinematic perturbations in section §6. A distance of 198 Mpc is assumed for 3C 120 
throughout the paper (Ho=50 km s~^ Mpc"^) which corresponds to a scale of ~ Ikpc/". 



2. Experimental Set-Up, Observations, and Data Reductions 

2.1. Integral Field Spectroscopy 

IPS optical observations with fibers are based on the idea of connecting the focal plane of 
the telescope with the spectrograph slit using a fiber bundle. In this way, when an extended 
object is observed, each fiber receives light coming from a particular region of the object. 
Each individual spectrum appears well separated on the detector and therefore, spatial and 
spectral information are collected simultaneously. The wavelength limitations essentially 
depend on the characteristics of the spectrograph itself. The spatial resolution depends 
on the fiber sizes and the prevailing seeing conditions during the observations. The spatial 
coverage of fiber systems is relatively small (see LeFevre et al. 2003, for a counter example), 
but they are very useful when studying small size objects, such us the circumnuclear region 
of nearby active galaxies (e.g. Garcfa-Lorenzo, Arribas, & Mediavilla 2001), blue compact 
dwarf galaxies (e.g. Cairos et al. 2002) or gravitational lenses (e.g. Motta et al. 2002; Wisotzki 
et al. 2003). 

IFS of 3C 120 was obtained the 26th of February 2003 at the Observatorio del Roque 
de los Muchachos (ORM) on the island of La Palma, Spain. The 4.2m William Herschel 
Telescope (WHT) was used in combination with the fiber system INTEGRAL (Arribas et al. 
1998, 1999; Mediavilla et al. 1998) and the WYFFOS spectrograph (Bingham et al. 1994). 
The observations were carried out under photometric conditions and an average seeing of 
l'/2. The standard bundle 2 of INTEGRAL was used during these observations. This bundle 
consists of 219 fibers, each one with a diameter of 0'/9 projected on the sky. A central 
rectangle is formed by 189 fibers covering a field-of-view of 16" x 12". 3, and the remaining 
30 fibers form a ring with a diameter of 90". Figure 1 illustrates the actual distribution of 
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the science fibers on the focal plane. 

The WYFFOS spectrograph was equipped with a 300 groove mm~^ grating centred on 
5500 A (spectral coverage: 3500-9000 A). A Tek6 CCD array of 1124 x 1124 pixels of 24 
//m size was used, giving a linear dispersion of about 3 A pixels"^. With this configuration, 
and pointing to the central region of 3C 120, three exposures of 1800 seconds each were 
taken. The data were reduced using IRAF standard routines (Tody 1986). Although the 
reduction of IFS data from fiber-based instruments docs not differ significantly from standard 
spectroscopic data reduction, we describe briefly in this section the reduction procedure. 

A master bias frame was built by averaging different bias frames taken along the night. 
This bias frame was then subtracted from the science frames. In observations with optical 
fibers, fiats-fields arc obtained illuminating the focal plane uniformly, and obtaining spectra 
(the so-called fiat-spectra). Thus, for a particular wavelength, the differences in response 
among fibers dominate the fiat-spectra shape. These differences are due to their distinct 
focal-ratio degradation, position at the entrance of the spectrograph, etc. In this way, flat- 
spectra are used to homogenize the response of all the flbers. Flat-spectra are also used 
to obtain the polynomial flts to deflne the flber path along the detector and extract the 
individual spectra from the whole image. Each spectrum appears well separated on the 
detector, with a width of approximately two pixels in the spatial direction (according to 
the core fiber image size of the fiber bundle used for these observations). The trace and 
extraction of individual spectra was performed using the standard routine APALL of IRAF. 
After this operation, frames of 1124 x 1124 pixels were reduced to 1040 x 219 pixel: each 
pixel in the spatial direction contains the spectrum of each particular fiber, a total of 219 
flbers. 

The wavelength calibration was done using the IDENTIFY and REIDENTIFY routines 
of IRAF. In order to carry out this procedure, we selected several isolated and well-distributed 
arc lines. Although it might be difficult to determine the actual uncertainty produced by 
the wavelength calibration, we have used the sky lines in our spectra to estimate the flnal 
wavelength errors, being smaller that 35.5 kms~^ and 29.5 kms~^ at blue and red wavelengths, 
respectively. 

The SP1045-I-378 flux standard star (Isaac Newton Group Database) was observed the 
same night and under similar conditions as for 3C 120. The standard star has been used to 
cahbrate a flux ratio by comparing with the standard-star flux tables of Stone (1977). We 
reduced the SP1045-I-378 frames in a similar way than the object ones, correcting them for 

differential atmospheric refraction (Filippenko 1982) using E3D (Sanchez 2004). Then, we 
extracted the observed spectrum of the star by co- adding the spectra of the central 37 fibers, 
which includes >99% of the total flux. Comparing this observed spectrum with the flux- 
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calibrated spectrum of SP1045+378 we derived a sensitivity curve that we used to calibrate 
the object frames. 

We combine the IDA tool (Garcia-Lorenzo, Acosta-Puhdo, & Megias-Fernandez 2002) 
and the Euro3D visualization package (Sanchez 2004) to analyze the data and generate two- 
dimensional maps of any spectral feature (intensity, velocity, width, etc). Maps recovered 
from spectra are images of 51x37 pixels with a scale O'.'3/pixel. While the spatial sampling 
of the used INTEGRAL configuration is 0'.'9 (that is, the fiber diameter of the bundle 2), 
the ccntroid of any peak in our maps can be measured with an accuracy of around 1/5 of 
the fiber diameter, that is ~ 0'.'2 (e.g. Mediavilla et al. 1998). 

2.2. HST imaging 

Wide-Field Planetary (WFPC) camera images of 3C 120 in different bands arc available 
on the archive of the Hubble Space Telescope (HST). We obtained these images to study 
the morphology of this object. Table 1 summarizes the properties of these images. The 
data set comprises three broadband images (F555W, F675W and F814W), which roughly 
corresponds to the standard V, R and /-bands, and a medium band image (F547M). These 
images basically sample the continuum emission in 3C 120, since the equivalent width of 
mostly all the emission lines is very small compared with the width of the bands. However, 
there is a non neglectible contamination, dominated by the [OHI] and H/3 emission lines in 
the F555W-band image (and Ha in the F657W-band image). The F547M-band image can 
be used as an estimation of the pure continuum emission, since the emission lines are at the 
edge of its transmission curve. The data set is deep enough to study the different structures 
present in this object (Soubeyran et al. 1989). 

3. Atleis of Spectra 

Figure 2 shows the nuclear spectrum of 3C 120 in the full wavelength range. This 
spectrum has been obtained by co adding the seven spectra closest to the continuum peak. 
This is almost equivalent to a hexagonal aperture of I'.'G in radius centred at the optical 
nucleus of 3C 120. We can easily recognize several emission lines and the characteristic 
broad component of permitted emission lines of Seyfert 1 galaxies. 

In the Appendix A, we present the individual spectra corresponding to each of the differ- 
ent observed positions (fibers) in selected spectral intervals that include the most important 
emission lines (spectra diagrams). 
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3.1. Emission-Line Profiles 

The profiles of the emission fines in tfie nuclear spectrum of 3C 120 sfiow a considerable 
blending due to the wings of the broad component of permitted lines (Fig. 2). Tfie H(5 
profile appears to be considerably broader tfian any otfier Balmer line, but tfiis is not tfie 
case of R5 as in previously found by other studies (Philfips & Osterbrock 1975; Baldwin 
et al. 1980). The HeIA5876 shows also a considerable broadening. At this step, we cannot 
discuss the origin of these features in terms of recombination models and reddening effects 
because the nuclear spectrum in figure 2 includes tfie contribution of tfie nucleus and tfie 
surrounding galaxy. In section §5 we separate botfi contributions and we will be in position 
to tackle tfiis discussion. Tfie low spectral resolution of tfie current IFS data prevents tfie 
detection of double peaks reported by Axon et al. (1989). However, tfie emission lines sfiow 
asymmetric profiles in some locations outside tfie central region. 

4. Data analysis and results 

4.1. Broad band distribution and morphological structures. 

Broad-band images of 3C 120 were recovered from our IFS data by coadding tfie flux in 
spectral ranges tfiat mimic tfie band-passes of tfie previously refereed HST images and using 
an interpolation routine (Sancfiez 2004; Garcia-Lorenzo, Acosta-Pulido, & Megias-Fernandez 
2002). The specific wavelengtfi ranges for eacfi band were 3900-4900 A(S-band), 5000-6000 
A(l/-band), 5600-5700 A(F'-band), 6100-7100 A(i?-band) and 7550-8550 A(/-band). AU 
images show a bright nucleus on top of a weak host galaxy. Figure 3A presents the in- 
tensity map corresponding to the y-band filter. Tfie intensity contours fiave an elliptical 
sfiape except for tfiose at around 4 arcsec from tfie central peak, sfiowing an elongation 
toward tfie west. Tfiere are remarkable similarities between tfie restored map (Fig. 3A) and 
the F555W-band image from the HST (Fig. 3B), despite of the differences in wavelength 
ranges and spatial sampling. The superior resolution of the HST images allows to directly 
detecting several structures. These structures have been previously reported subtracting a 
galaxy template to broad and narrow band ground base images (e.g. Baldwin et al. 1980; 
Perez-Fournon et al. 1986; Hua 1988; Soubeyran et al. 1989). We have labeled, using the 
nomenclature introduced by Soubeyran et al. (1989), the different detected structures (figure 
3B). Tfie HST/F555W-band image allows to resolve tfiem, sfiowing a ricfier level of struc- 
tures. In particular, structures B, and C are composed of a smootfied low surface brigfitness 
component and more luminous clumpy substructures. Structure A seems to be more colli- 
mated, as previously noticed by ground based imaging, wfiich explains wfiy it was initially 
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confused with the optical counterpart of the radio jet (e.g. Perez- Fournon et al. 1986). It 
shows also a low surface brightness component and four clear clumps. We overplot in figure 
3B the radio-map at 4885 MHz (Walker 1997), showing that the structure A is most likely 
related to the structure B rather than to the radio-jet. Indeed, the HST images suggest 
that both structures belong to a sequence of clumpy knots which are physically connected. 
Dust lines and shell structures not detected from the ground are seen in the inner regions. 
Two shell arcs at ~ 1" north and south of the central peak are clearly detected in the HST 
images. Although previous HST images show extensions to the north-west and south-east 
in the inner region of 3C 120 (Zirbel & Baum 1998), to our knowledge this is the first time 
that the shell structure is clearly visible. The north and south shells are named S;v and S5, 
respectively, hereafter. The S5 seems to be connected to structure A by a faint tail. 

A 2D modehng of the broad-band images was performed using GALFIT (Peng et al. 
2002) to obtain a clear picture of the morphology of 3C 120. This program has been ex- 
tensively tested in the image decomposition of QSO/hosts (Sanchez et al. 2004b). The 2D 
model comprises a narrow Gaussian function (to model the nucleus) and a galaxy template 
(to model the galaxy) both convolved with a PSF. We performed the fit twice, first using 
a Sersic law (Sersic 1968) to characterize the galaxy and then using a de Vaucouleurs law 
(de Vaucouleurs 1948). This method allows us to determine the morphological type of the 
host galaxy, based on the obtained Sersic index, and to get a good determination of the 
galaxy flux, based on the modehng the de Vaucouleurs law (Sanchez et al. 2004b). The PSF 
was created using a field star for the analysis of the HST images, and a cahbration star for 
the analysis of the broad-band images recovered from the IFS data. 

We applied first the fitting technique to the HST images because of their better spatial 
resolution. The nucleus was saturated in the HST/F814W-band, which prevents us to do a 
two-model fitting, and limits the reliability of any morphological classification, which strongly 
depends on the shape of the profile in the inner regions. To derive a rough estimation of 
the host magnitude we masked the nucleus, and fixed the scale to the average of the values 
derived in the other available HST/bands. Further checks, explained below, demonstrated 
that this approach was valid. Table 2 summarizes the results from this analysis. For each 
band it shows the derived Sersic index, the nucleus and host magnitudes and the effective 
radius of the host galaxy. The derived Sersic indices, all near or larger than the nominal 
value of 4 for an early-type galaxy, confirm the morphological classification of the host galaxy 
of 3C 120 as a bulge dominated galaxy. 

We subtracted the object template (galaxy+nucleus) derived by the 2D modeling from 
the original images, obtaining a residual image for each band. These residual images were 
used to study the properties of the different structures, once decontaminated from the smooth 
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component. As an example, the residual image of the F-band and the HST/F555W-band 

are shown in Figure 3C and 3D, respectively. The already quoted structures are now clearly 
identified. The S5 shell coincides with an extended emission line region (EELR) previously 
detected in this object (Hua 1988; Soubeyran et al. 1989; Sanchez et al. 2004a) and labeled 
as E2 (Soubeyran et al. 1989). 

4.2. Colors and gaseous distributions in 3C 120 

As already quoted all the broad-band images recovered from the IPS data show a similar 
morphology than the V-band image shown in figure 3A. Combining these recovered broad- 
band maps, we obtained color maps of 3C 120. Figure 4A shows the V — I color map derived 
from the IPS data. It shows a blue circumnuclcar region elongated to the southwest. The 
elongation towards the west seen in the broad band images has also a counterpart in the 
color maps showing bluer colors at a region ~ 4'.'5 west of the nucleus than its surroundings. 
Figure 4B presents the V — I color map derived from the F555W and F814W-band HST 
images. Despite of their different spatial resolution, figures 4 A and 4B present remarkable 
similarities. The different structures quoted in section §4.1 show bluer V — I colors than the 
average color of the host galaxy, as already noticed (e.g. Praix-Burnet et al. 1991). Figure 
4B also shows evidence for two weaker shells farther from the central peak than and S^ 
but also at the north and south. However, these "secondary" shells are fainter, being at the 
detection limit. 

Both the y-band and the F555W-band images are contaminated by the emission from 
[0111] and H/3. A rough estimation of how strong that contamination is can be obtained 
by subtracting a scaled continuum image, clean of those contaminations. Figure 4C shows 
the residual image of the INTEGRAL l/-band image after the subtraction of an adequate 
continuum (the V'-hsaid image). A similar estimation was done for the HST images, sub- 
tracting the continuum dominated F547M-band image from the F555W-band image (Figure 
4D). This image was smoothed using an 11 x 11 pixels median kernel to increase the signal-to- 
noise ratio. The resulting maps are a rough estimation of the distribution of the [0III]-|-H/3 
emission. The position of two of the EELRs detected in this object (Ei and E2, hereafter) is 
indicated (Soubeyran et al. 1989; Sanchez et al. 2004a). The residual images of the contin- 
uum, after subtracting a model template (Fig. 3C and 3D), are overplotted for comparison 
purposes. The Sjv and S5 shells show a strong gaseous emission in the HST images, which 
in the S^ extends towards the E2 region. 

Figure 5 shows the V — R colors of the host galaxy and the different structures as a 
function of the R — I colors. As already noticed in the color image, the structures have 
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bluer V — R colors than the host galaxy. This indicates most probably the presence of 
somewhat younger stellar populations associated with them. However, they show a wider 
range of R — I colors than that of ^ — i? colors. Indeed, the average R — I color of the 
host galaxy is within the range of i? — J colors of the different structures. The color-to-color 
distribution derived from synthetic models has been included in the figure for comparisons. 
The discontinuous line shows the colors of single stellar populations calculated using the 
Bruzual & Chariot (2003) models. We assumed a solar metahcity and a Chabrier (2003) 
IMF. The labels indicate the logarithm of the stellar population ages in Gyrs. In general 
terms, the color-to-color distribution does not match with that simple model, apart from 
the case of the Sjv and S5 shells. For the remaining structures, the R — I colors correspond 
on average to a young stellar population of ~6xl0^ Gyr, but the V — R colors correspond 
to older stellar populations. This may indicate that the real populations are a composite of 
populations. The dotted line shows the colors of different galaxy types (E, Sab, Sbc, Scd and 
Irr) (Fukugita et al. 1995), which mainly correspond to different mixes of stellar populations. 
The R — I color of the host galaxy mainly corresponds to a Sab-Scd galaxy. On average, the 
structures have V — R colors that correspond to the same kind of galaxies (Sab-Scd), but 
their R — I colors expand over all the range of possible colors. 

So far we have ignored the effects of dust in the color-to-color distribution shown in 
Figure 5. However, the dust content in 3C 120 is rather high; an average value of Ay ~ 4 
mags has been estimated (Sanchez et al. 2004a). Dust extinction redden both V — R and 
R — I colors in an almost similar way (Fitzpatrick 1999). The effects of dust in the color- 
to-color distribution is illustrated with an arrow in Fig. 5. It is clear that the dust is not 
homogeneously distributed in the galaxy: e.g., typical dust lanes are seen in figure 3B. A 
combined effect of different stellar populations and a non-homogeneous distribution of dust 
could explain the observed color-to-color distribution. In that case, the Sjv and S5 structures 
would be dust-free areas dominated by a single stellar population. 

However, results from broadband colors should be taken with care. The contamination 
of broadband filters by strong emission lines can drastically affect the morphology of color 
maps. In the case of 3C 120 the problem is even worse because of the contribution and 
contamination from the broad component of the Balmer emission lines in the central regions. 

4.3. Line-Intensity Maps and Ionization Structure 

We performed a line profile fitting and de-blending in order to study the integrated high 
and low ionization gas distribution, their physical properties and kinematics. We fit a single 
Gaussian to any of the bright emission lines in the spectra. For the Balmer lines, we included 
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a second broader Gaussian to fit the wide line from the broad hne region (BLR). In spite of 
our poor spectral resolution, we found that several spectra have evidence of substructures, 
showing asymmetric profiles with blue or red wings indicating the presence of several gaseous 
systems in 3C 120. Indeed, previous authors (Axon ct al. 1989) reported the existence of 
several components in the emission lines of 3C 120, but to apply a line-profile decomposition 
to separate the different gaseous components would be unrealistic because of the low spectral 
resolution of the current IFS data. The broad component of the Balmer lines is confined to 
the nucleus discarding an extended broad emission line region. 

Figure 6 shows the intensity maps of [OIII]A5007 (figure 6A) and Ha (figure 6B) after 
the deblending. Intensity contours of ionized gas are clearly elongated toward the E2 struc- 
ture. Emission line maps present an elongation to the west, forming a clear secondary peak 
at around 5 arscec from the nucleus in the [OIII] intensity map. This secondary peak corre- 
sponds to the El structure previously quoted (section §4.2). The ratio of the Ha intensity 
of the nucleus and Ei is much smaller than that of [OIII]A5007, suggesting a high ionization 
nature of the latter. The nuclear region shows a low ionization degree in the [OIII]A5007/H/3 
map surrounded by a ring of high ionization (figure 6C), with larger values at the south of 
the nucleus, coincident with the location of E2. A high ionisation region expands from the 
ring to the west, increasing the ionisation degree along this direction. The [NII]A6584/Hq; 
line ratio map (figure 6D) also shows a ring structure, surrounding the nuclear region. The 
El is on the path of the 3C 120 radio-jet, as well as the continuum dominated structure A. 

The high [0III]/H/3 ratio of Ei (figure 6C) discards a relation between A and Ei, as- 
suming that A is a. star forming region in a spiral arm or a tidal tail, and points, more 
likely, to a direct connection between Ei and the radio-jet (Sanchez et al. 2004b). While the 
nucleus and the A region present [OIII]A5007/H/3 and [NII]A6584/Hq; ratios at the limit of 
HII regions in a diagnostic diagram (Veilleux & Osterbrock 1987), those ratios of Ei and its 
west surroundings are placed on the high-ionization region . 

The distribution of the Ha/H/? hne ratio (figure 6E) shows a dust lane crossing 3C 120 
along the southeast-northwest direction. The high Balmer decrement at the west points to 
a high dust obscuration around Ei. However, we cannot rule out that the high Hq;/H/3 ratio 
may be related to the low signal to noise of the narrow component of H/3 in the circumnuclear 
region of the galaxy. The electronic density derived from the [SII]A6716/A6730 presents a 
patchy structure, showing an enlargement in the location of Ei. In the next section, the 
ionization conditions of the different structures in the observed region of 3C 120 are studied 
in detail. 



5. Decoupled spectra of the different components in 3C 120 



In the previous section, we used the traditional method (Gaussian fitting) to deblend 
the different gaseous components. IFS optical observations record spatial and spectral in- 
formation simultaneously. Taking advantage of this fact, we have developed a technique 
(3D-modeling hereafter) to disentangle the spectra of the main components of 3C 120 (nu- 
cleus+host) , and to obtain cleaned spectra of the different structures described above. This 
technique has been successfully apphed recently (Sanchez et al. 2004a). A brief summary of 
the technique is described in the Appendix B, and will be explained in detail in a separate 
article (Sanchez et al., in prep.). In summary, the technique provides us with a spectrum 
of each of the main components of the object plus a residual data cube that can be used 
to derive cleaned spectra of the different structures in the object and/or analyze their mor- 
phological and kinematical properties. This is an extension of the 2D modeling of images 
applied in section §4.1 to each of the monochromatic images of the data cube derived from 
IFS optical observations. We apphed the 3D modeling to the current data and we derived 
the 3C 120 nucleus spectrum, the mean host galaxy spectrum and a residual data-cube of 
spectra. 



5.1. Morphological structures from the residual data cube of spectra 

To compare the results from the proposed technique to those of the traditional method, 
figure 7A shows the greyscale of the residual from the 2D modehng of the F555W/HST-band 
image obtained in section §4.1 (figure 3D). We over plotted the contours of a narrow-band 
image centred on the continuum adjacent to the [OIII] emission fine (5204-5246 A) extracted 

from the residiial data cube of spectra obtained after the 3D modeling. The main continuum 
dominated structures (A, B and C) are detected ~1" from the nucleus. As expected, the 
shell structure detected at ~1" north and south of the nucleus in the HST images is not 
detected with the IFS data. A combination of the superior resolution of the HST images, 
the effects of the seeing and the sampling of our IFS data can explain it. In figure 7B we 
present the greyscale of the F555W— F547M image (obtained in section §4.2, figure 4D) and 
the contours of a narrow-band image centred on the [OIII]A5007 emission line (5170-5200A) 
extracted from the residual data cube. As we quoted above (section §4.2), the F555W-F547M 
image is a rough estimation of the [0III]-H/3 emission map. Although both images are in good 
agreement, integral field spectroscopic data are more adequate to detect extended emission 
line regions (EELRs) that are diluted in broadband filter images, like the F555W-band one. 
It is possible to identify the Ei and E2 structures, and a third region (E3 hereafter) at the 
north-west of the F555W-F547M image. Soubeyran et al. (1989) already detected these 
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emission line regions, using narrow-band imaging centred on the [OIII] emission line spectral 
region. Ei and E2 are completely coincident with their reported positions. The peak of E3 
found by Soubeyran et al. (1989) is not included within the field of view of INTEGRAL. 
Wc arc confident of our detection, and most probably we are seeing a tail previously not 
detected of E3 towards the nucleus. 

5.2. Spectra of the nucleus and the host galeixy 

Figure 8 and 9 show the spectra of the different deblended components of 3C 120. Figure 
8A shows the spectrum of the nucleus. Its continuum follows roughly a power-law with a 
slope of a ~ —0.5 {Fx oc A"). It contains broad-emission lines, without any trace of them in 
the spectrum of the host galaxy (figure 8B), which indicates that the decoupling technique 
has worked properly. The intensities of the narrow emission lines is weaker in the nucleus 
than in the host galaxy spectrum. Indeed, many emission lines are clearly detected in the 
host galaxy spectrum but not in the nucleus. 

To derive the properties of the emission lines in the spectrum of the nucleus, we fit 
Gaussian functions to the different emission lines. The adjacent continuum was modelled 
with a low-degree polynomial function. To increase the accuracy of our model, different 
nearby emission lines were fitted together, defining line systems with the same systemic 
velocity and full- width-half-maximum (FWHM). After several attempts, the best (and sim- 
plest) modeling was obtained including two decoupled broad systems, with ~8796 km s~^ 
and ~2188 km s^^, and a narrow system, with ~700 km s~^. Recently, Ogle et al. (2004) 
also reported two broad and one narrow component for the H-Balmer lines in the nucleus 
of this galaxy with a similar velocity dispersions than those derived here. Table 3 lists the 
result from the fit, including the wavelength, the fiux and their la errors for each detected 
line. The uncertainties in the determination of the flux of the narrow emission lines affected 
by fine blending are high as expected. The Ha/H/? fine ratios are in the range of ~3.0-3.6 
for the different line systems, in agreement to the nominal values for case B recombination 
(Osterbrock 1989). Therefore both the narrow and broad emission line regions are not sig- 
niflcantly affected by dust, as expected for a type 1 AGN nucleus. We estimated an effective 
temperature of ~2.3xl0-'^ K using the ([OIII]A5007+[OIII]A4959)/[OIII]A4363 line ratio and 
the relation between this ratio and temperature (Osterbrock 1989) typical of this kind of 
AGN. 

Figure 8B shows the average spectrum of the host galaxy together with the spectra 
of different synthetic models consisting on single stellar populations of different ages. We 
used the Bruzual & Chariot (2003) models, assuming a solar metallicity and a Chabrier 
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(2003) initial mass function (IMF). The effects of dust were included in the model using the 
extinction curve of Fitzpatrick (1999) for a dust absorption of Ay ~4 mags (the average value 
in the host galaxy, as we will show below). Two extreme cases were plotted, representing an 
old (~16Gyr, red line) and a young stellar population (~0.01Gyr, blue line). These extreme 
clear over simplification, but they can be used just for qualitative comparisons. A 
single stellar population cannot describe the observed spectrum: the optical slope between 
A ~5000 A and A ~6500 A matches well with that of an old stellar population. However, the 
slope at larger wavelengths is flatter, more similar to that of a young stellar population. We 
already noticed it (section §4.2) when we analyzed the broadband colors of the host galaxy 
(Figure 5). Only a mix of different stellar populations can explain the observed spectrum. 

We included in figure 8B the spectrum of a fifty-to-fifty mix of both components (orange 
line). Even this simple model can describe quite well the spectrum of the host galaxy for 
A >4500A. At shorten wavelengths there is an increase of the fiux, which does not match with 
the mix model (dashed line). It is well known that powerful radio galaxies show a significant 
blue-UV excess (Lilly & Longair 1984) most probably due to scattered light from the nucleus 
than to young stellar populations (McCarthy et al. 1987). This scattered light can be well 
described by a power-law. Adding a power-law to the mix model we can describe the average 
host galaxy spectrum for any wavelength (orange solid line). Therefore, three components 
are needed to explain qualitatively the spectrum of the host galaxy: (a) an underlying 
old stellar population, (b) a young stellar population, and (c) a nebular continuum, most 
probably scattered light from the nucleus. We did not perform a proper fitting to determine 
the best fractions of the different components because it would not change the qualitative 
statement, and, in any case, the specific fractions will depend on the assumed synthetic 
models for the young and old populations and we chose them arbitrary. The contribution 
of each component to the total flux are ~5% (a), ~55% (b) and ~40% (c) at 4500 A and 
~35%(a), ~64% (b) and ~2% (c) at 6000 A. Similar results have been found in the study 
of other powerful radio galaxies (e.g. Tadhuntcr ct al. 1996). 

Table 3 lists the results derived from fltting a single Gaussian to the emission lines in 
the host galaxy spectrum. As we quoted above, only narrow emission lines corresponding 
to a dispersion velocity of FWHM~700 km s~^ are detected. These emission lines are 
considerably stronger than the narrow-emission line region in the spectrum of the nucleus, 
indicating the presence of an extended narrow-emission line region which extends through 
most of the field of view. Similar results have been found in spectroscopic studies of radio- 
quite type 1 AGNs (e.g. Jahnke 2002; Jahnke et al. 2004). The Ha/H/? fine ratio is ~10, 
indicating a dust absorption of Ay ~4 mags when comparing with the nominal case B 
recombination value (Osterbrock 1989) and using the Fitzpatrick (1999) extinction curve. 
In average the host galaxy contains a rather high dust content, but the dust distribution is 
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not uniform as we can see in figure 6E. We estimated an effective temperature of ~14000 
K using the ([OIII]A5007+[OIII]A4959)/[OIII]A4363 fine ratio and an electron density of 
ne~10 cm^^ using the [SII]A6716/6731 Unc ratio. To derived these quantities, wc used the 
relations between these line ratios and the measured parameters (Osterbrock 1989). The 
fine ratios log([Olll]A5007/H/3)~1.2 and log([NII]A6583/Ha)=-0.4 indicate most probably 
a direct photo ionization by the UV flux from the AGN, consistent with the result shown in 
Section 54.3. 



5.3. Spectra of the EELRs: Origin of the ionization 

The integrated individual spectra of the structures Ei, E2, E3 described in section 
§5.1 were extracted from the residual data cube. Figure 9A presents the spectra of each 
EELRs (El, E2, E3), which are remarkable similar. They have a clear gaseous nature, 
without appreciable continuum emission. This result confirms that the continuum dominated 
structure A is not related to Ei as we pointed out in section §4.2. A visual inspection of the 
relative strength of the different lines indicates that the dominant ionization source is the 
AGN or shocks, rather than a star formation process. This was already noted in section §4.3 
and figure 6C for Ei, and E2. Table 4 lists the result of the Gaussian modeling of the emission 
lines in these spectra. Different line ratios and parameters derived from them are listed in 
Table 5. We included the [Olll]A5007/H/3, [N11]A6583/Ha, [S11]A6716/[S11]A6731 line ratios, 
and the derived dust absorption {Ay), electron density (ue) and effective temperature (Tg//). 
Those values of the nucleus and the host galaxy spectra were also included for comparison 
purposes. As we already mentioned, the ionization conditions are similar in the different 
EELRs, being also similar to those of the nucleus and the host. Figure 10 shows the classical 
diagnostic diagram of the [OIII]A5007/H/3 line ratio as a function of [NII]A6583/Hq; for the 
different components of 3C 120, including the division between AGNs and star forming 
regions (Veilleux & Osterbrock 1987). The major difference is found in the position of the 
nucleus in the diagram, which seems to lie in the location of star forming regions. This is 
most probably due to the large uncertainties in the [N11]A6583/Ha! line ratio derived for the 
spectrum of the nucleus, largely influenced by the deblending process of the narrow and broad 
emission fines, as we quoted above. The derived effective temperature (T~ 1.4 x 10^ K) for 
the host galaxy spectrum does not exclude the hot stars as a possible origin for the ionization 
but the high [SII]/Hq;, [NII]/Hq;, and [0III]/H/3 ratios indicate that the average ionization 
for the host is the AGN or/and a shock process, i.e., the observed spectra cannot be found 
in HII regions. The [N11]A6583/Hq! line ratio of the E3 region is lower than that value for the 
rest of the structures. E3 is the faintest detected EELR and it was not completely covered 
by our field-of-view, according to the [OIII] maps shown by Soubeyran et al. (1989). Its 
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spectrum is noisier and with some clear defects at specific wavelengths produced by the 
modeling technique. A detailed inspection of the spectrum shown in Fig. 9 shows that the 
[NII]A6583 line is distorted by one of these defects, affecting the [N1I]A6583/Hq; line ratio, 
which error is clearly larger than the formal error plotted in Fig. 10. Despite these caveats, it 
is clear that star formation processes do not dominate the ionization of the different EELRs. 

The dust content in the EELRs is lower than the mean obscuration of the host. In 
particular, for the E2 and E3 regions the dust absorption is almost half of the average in the 
host galaxy. On the other hand, the electron density in those clouds is ~10 times lower than 
the average. However, this density is ~10 times higher in the Ei region, i.e. ~100 times 
larger than in the other clouds. This indicates most probably a different origin of those 
clouds. In Sanchez et al. (2004a) we discussed in detail the nature of the Ei emission line 
region. It is most probably associated with the radio jet that crosses this cloud, compressing 
it due to its lateral expansion, and splitting it in two different kinematics regions (Axon 
et al. 1989; Sanchez et al. 2004a). The compression is reflected in the increase of the electron 
density. The nature of the ionization is not clear in this cloud, since a post-shock zone 
can also give rise to the observed line ratios. The similarities between the line ratios for 
the different EELRs, and between them and the host galaxy ones, may indicate a similar 
origin for the ionization. That is, direct photo ionization from the UV-field of the AGN. In 
that case the effect of the jet over the Ei would be reduced to a split of the cloud and a 
compression that give rise to a density enhancement. 

5.4. Spectral energy distribution of the continuum dominated condensations 

Figure 9B shows the spectral energy distribution (SED) of the continuum dominated 
structures A, B and C detected in 3C 120 within the field-of-view of our IFS data. The SEDs 
were obtained by an average of the integrated spectra of the different structures extracted 
from the residual data cube over spectral ranges of SOOA width. The low intensity of the 
continuum dominated structures and the subsequent low signal-to-noise prevents us of using 
directly the extracted spectra for this analysis. The SEDs were cut at 7500A since at 
larger wavelengths the imperfect subtraction of the sky-lines and the noise enhancement 
strongly affect the reliability of the derived SED. Synthetic spectra for three different stellar 
populations of different ages were included for comparison purposes. These spectra are 
similar to those shown in the figure 8B (described in §5.2), but without including the effects 
of dust. The SEDs are almost flat in the plotted ranges, being roughly consistent with a mix 
of a young stellar population with an underlying old stellar population. This result agrees 
with the results based on the broadband colors of the structures, discussed in §4.2. Therefore, 
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the continuum structures are experiencing a decrease of the dust content, rather than an 
increase of the star formation. This seems to be particularly valid for the A structure, which 
lies in a minimum of the dust content, derived from the Ha/H/? line ratio (See Fig. 6E). 
Unfortunately, we cannot check it for the B and C structures, because the gaseous emission 
in those regions is less luminous and therefore the Hq;/H/3 ratio is uncertain. 

So far only starlight was considered to explain the emission found in the different con- 
densations. However their SEDs could be also explained considering other components, 
like scattered-light from the nucleus and/or synchrotron emission associated with the radio- 
jet. Hjorth et al. (1995) detected polarization in the condensation A, which direction and 
magnitude were consistent with those found by Walker et al. (1987). This may indicate a 
connection between that condensation and the radio jet. Indeed, it may indicate that con- 
densation A consists, at least partially, of optical synchrotron emission. However, as it was 
shown in §4.1 (and noticed by Hjorth et al. 1995), A does not follow the radio-jet into the 
core, being more likely associated with B and C. Furthermore, the measured polarization 
may be due to unsubtracted scattered-light from the nucleus (Hjorth et al. 1995), that we 
detected in the average spectrum of the host (Fig. 8 and §5.2). Indeed, the SEDs of the 
three different condensations do not differ significantly (Fig. 9B), and, in particular, they 
show almost the same V — R colors (Fig. 5), which indicates, most probably, a similar origin 
for all of them. Since most of the scattered-light has been removed and the synchrotron 
radiation could also contaminate A, these condensations are most probably dominated by 
starlight. 



6. 3C 120 Velocity Field 

The central wavelengths of the Gaussians fitted to the individual spectra (section §4.3) 
give us the radial velocity associated to the ionized gas at each position. Interpolating the 
individual emission line centroids, we obtain the velocity field of different emission lines. 
Uncertainties due to wavelength calibration (< 35 km s~^) and those related to the fitting 
process (~ 15 km s~^) are small enough to be irrelevant for the following discussion. 

Figure 11 shows the velocity field of the ionized gas derived from [OIII]A5007 and Ha 
fines. While Ha traces the kinematics of low ionisation gas that describe the general pattern 
of the galaxy, [OIII] draws the signatures of high ionization gas, characterizing the most 
perturbed regions. Previous kinematics studies of 3C 120 found a rather chaotic velocity 
field (Baldwin et al. 1980) and some slight evidences of a co-existing rotating system (Baldwin 
et al. 1980; Moles et al. 1988). Despite of the different spatial resolution and coverage of the 
IFS data in this paper, our results are in agreement with those of Baldwin et al. (1980) and 
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Moles et al. (1988) pointing towards highly distorted kinematics in this object. 

The unclear morphology and highly distorted gas kinematics of 3C 120 makes a simple 
interpretation of the velocity structure difficult. The ionized gas velocity field presents a 
general regular pattern, with larger velocities at the north-west and lower at the south-east 
that may be consistent with a disk rotating around an axis along the north-east direction 
(figure llA and IIB). Several kinematical perturbations can be identify in the velocity maps 
as well as a regular rotation. 

A velocity gradient is located at ~5" west of the nucleus, aligned with the radio-jet 

(figure 3B). at the position of Ei (figure 11 A). Using a better spectral resolution, Axon 
et al. (1989) show that there are two different kinematics components rather than a velocity 
gradient. The north component is receding while the south component is approaching. 
The kinematical perturbation is most probably due to the lateral expansion of the radio-jet 
(Axon et al. 1989): as already quoted, the interaction of the radio-jet and the intergalactic 
gas produces an enhancement of the gas density (section §4.3), and perturbs the kinematics 
(Sanchez 2004). 

The velocity field derived from [OIII] is much more distorted to the east (farther than 4 
arcsec from the nucleus) than that one derived from Ha. Although the signal to noise in the 
IPS spectra is low and uncertainties are larger in that region, differences in the kinematics 
may be explained by the different origin of these lines. The comparison of the continuum 
dominated structures and the velocity field derived from the Ha shows that A and B coincide 
with a velocity distortion located at the west of the nucleus (figure IIB). This kinematical 
feature is similar to those found in the velocity field of spiral galaxies because of the arms 
(e.g. Knapen et al. 2004). This gives support to the idea of identifing the A and B regions as 
structures in an spiral arm in 3C 120 in spite of its poor gas content and continuum emission 
dominated by starlight. We will discuss in detail the origin of this arm-like structure in 
section §7. 

The comparison of the ionized gas structures E2, and E3 (figure 7B) and the velocity 
field gives clues to the origin of the EELRs (figure IIC). The morphology of E3 follows 
remarkably well an east-west gradient in the velocity map. E2 is in a close to constant velocity 
region. Both clouds are in regions which present strong kinematical perturbations from the 
canonical rotation. These perturbations are more prominent in the velocity structure of the 
high ionization gas than that of the low ionization gas. 

Fitting a single Gaussian to the emission lines in the residual data cube of spectra 
(section §5), we have derived the velocity behavior of the Ei, E2, and E3 residual structures. 
Figure IID shows the velocity map derived from the [OIII] hues in the residual spectra. 
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We will refer to this map as the velocity map of residuals hereafter. The kinematics of 
the structures trace remarkable well the perturbations in the velocity field, although their 
kinematical features are smoother in the latter. The smoothing is a clear consequence of 
the blending of the host galaxy and the structure spectra. The velocity map of the residuals 
represents much better the kinematics of the gas clouds in Ei, E2, and E3. We found a 
clear gradient of velocities in the region corresponding to the three structures pointing to 
inflows/outflows of gas at different inchnation angles and outside of the galactic disk. While 
regions Ei and E2 present a north-south velocity gradient, E3 shows an east west gradient 
in good agreement with the [OIII] velocity map (figure llA). 

Despite of the bulge dominated morphology found in the inner regions of 3C 120, the 
velocity field derived from Ha indicates the presence of a rotating disk along the north- 
west/south-east (figure IIB). In that case, the receding velocities at the northwest suggest 
that southeast is the face closer to the observer. The velocity gradient also traces the hne 
of zero velocities, but it is not clearly defined in the derived velocity field because of the 
distortions in the center introduced by the Syl nucleus. With a larger spatial coverage, 
Baldwin et al. (1980) determined a PA=72° ± 15° for the minor kinematics axis. This angle 
agrees with the semi-minor axis determined from the 2D modeling of the galaxy (section 
§4.1). A visual inspection of the velocity fields in figures llA and IIB indicates that this 
PA may be a good estimation of the minor kinematics axis, when considering the distortion 
at the west. However, this distortion is most likely produced by the interaction of the radio- 
jet with the El cloud (Sanchez 2004), and most probably is not related with the rotating 
disk. The isovelocity contours of the velocity field derived from Ha and the position of the 
kinematical center (~ 0.8" north-west of the nucleus) suggest that the minor kinematical 
axis is most likely along PA~ 50°, which is not far from the photometric determinations and 
within their uncertainties. Therefore, we will consider hereafter a PA=— 40° for the major 
kinematics axis of 3C 120. 

Figure 12 shows the rotation curve of 3C 120 derived from the velocity field of Ha, 
assuming a PA of —40° for the major kinematics axis. This curve is clearly distorted at the 
south-east due to the kinematical perturbations associated with E2. However, the north- 
west portion of the curve is remarkable similar to those of galactic rotating disks (Binney 
& Merrifield 1998). The effects of the perturbations in the south-east portion of the curve 
must be removed prior to model the 2D distribution of the rotating component of the galaxy. 
For doing so, it was assumed that the approaching portion of the rotation curve (south-east) 
follows a symmetrical counterpart of the receding portion (north-west). This hypothetical 
rotation curve was included in Figure 12. This curve was then used to derive a template 
of the velocity field by applying a simple rotational model (Mihalas & Binney 1981), for 
different inclinations varying from ~5° to ~85° . Subtracting those templates from the Ha 
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velocity field and minimizing the differences at the northwest region (where the rotation curve 
really corresponds to the 3C 120 kinematical behavior), we estimated an inclination angle of 
~ 40° for the disk component of 3C 120. The velocity field template corresponding to that 
inclination angle is shown in figure HE. This estimation of the inclination is in agreement 
with the determination from the external isophote (Moles et al. 1988). As quoted before, the 
internal isophotes are less elongated and distorted at the west than the external ones. This 
suggests a different inclination angle, smaller at the inner than at the outer regions. Indeed, 
the ellipticities derived from the 2D modehng (section §4.1) indicate that the inchnation 
may range from ~ 20° to ~ 43° in the center, in agreement with the kinematical estimation. 
The differences in inclination from the center to the outer regions, although small, suggest 
a slightly warped disk component in 3C 120. 

Figure IIF shows the residual of the [OIII] velocity map once subtracted the template 
shown in Fig. HE (residual velocities hereafter). This residual map shows an almost flat 
area of 0±20 km s~^at the north-west region, within the uncertainties of the velocity de- 
termination. At the west, coincident with the path of the radio-jet and the location of Ei, 
there is a north-south gradient of ~ 148 km/s, in agreement with the results by Axon et al. 
(1989). At the location of E2 there is a rather chaotic residual velocity structure, with veloc- 
ities ranging from ~270 to 460 km/s. A slight east-west gradient of ~ 32 km/s is found at 
the location of E3. The residual velocities (Fig. IIF) and the velocity map of the residuals 
(Fig. IID) present a rather good quahtative agreement in spite of the strong conceptual and 
practical differences of the two methods used to derived them. This supports the idea of the 
presence of a rotational component in 3C 120. 

Due to the low spectral resolution of our data no attempt was done to analyze the 
velocity dispersion maps. 

7. discussion 

In previous sections wc presented several aspects of 3C 120 that describe a puzzle en- 
vironment. 3C 120 is morphologically speaking a bulge dominated galaxy, which contains a 
rotating stellar disk and several continuum- dominated structures and EELRs. We analyzed 
in detail those structures in order to determine their nature. The Ei is caused by the inter- 
action of the radio-jet with the intergalactic medium, as already discussed in Sanchez et al. 
(2004a). 

The newly reported shells in the central region of 3C 120 is an addition to the complex 
picture of this object. Shells are common structures in elliptical and SO galaxies, although 
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a few number of spirals also show the presence of this kind of features (Mahn & Carter 
1983; Schweizer & Seitzer 1988). Although models considering internal shock waves were 
proposed to account for the presence of shells in galaxies (see e.g. Williams & Christiansen 
1985) the most accepted idea is that they are generated by merging processes (see e.g. 
Hernquist & Spergel 1992, and reference therein). Numerical simulations indicate that the 
number and sharpness of shells in merger remnants depend on the mass ratio and the absence 
(or presence) of a central buldge in the progenitors (Gonzalez-Garcia & Balcells 2004). The 
radial distribution of shells depends on the potential of the host galaxy and their morphology 
can appear aligned or randomly distributed around the galaxy (Prieur 1990). 3C 120 shows 
two well defined shells, labelled as Sjv and S5, and some other shell-like faint features which 
are aligned along an axis of PA~ —25°. When shells are well aligned along a certain axis, 
this axis is always close to the major axis of the galaxy (Prieur 1990), which is out case. 
Shells can be found in a wide range of radii from the nucleus, from only a few to hundreds of 
Kpc, such us in the case of NGC 3923 (Prieur 1988). The innest shell detected in 3C 120 is 
at around 1 arcsec north (~ 1 kpc) of the nucleus. Shells found close to the nucleus implies 
that a dissipative process has played an important role in the formation of the shells (Prieur 
1990). Therefore, the presence of the detected circumnuclear shells (section §4.1) strongly 
suggests the idea of considering 3C 120 as a late stage merger. 

The arm-like clumpy structures A and B (section §4.1) show kinematics that resembles 
those of a spiral arm as explained in section §6. However, spiral arms are active star forming 
regions, and no star formation was detected in either of the continuum-dominated structures. 
In that case, these structures could be also the remnant of a merging process. Indeed, the 

regions B, A and the Sg seem to be connected in our residual and color maps (Fig 3C and 
4B), and S^ is associated with the E2 EELR. A kinematics feature (Fig. IIB) is coincident 
with these structures. Numerical simulations indicate that these kind of structures are found 
as remnants of merging processes (e.g. Howard et al. 1993; Heyl, Hernquist, & Spcrgcl 1996; 
Mihos & Hernquist 1996) being difficult to be explained by any other mechanism. Assuming 
that scenario, the presence of an over density of gas in the south shell may indicate an 
infiow, where the warmer gas has dropped faster into the inner regions. This picture is in 
good agreement with the kinematics features at the south shell (E2 clouds kinematics) and 
its residuals. The regular velocity in this region suggests that this inflow may occur in a 
plane outside the galactic disk and close to the line of sight. Although the kinematics of 
the merger remnants strongly depends on the viewpoint, the velocity and velocity dispersion 
measured for E2 is in good agreement with merger simulations (Heyl, Hernquist, & Spergel 
1996). Theoretical models indicate that inflows channeling gas from the outer to the nuclear 
regions may appear in the late stages of a merging process (e.g. Mihos & Hernquist 1996). 
These inflows have been already observed in different merging galaxies (e.g. Arribas & Colina 
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2002, 2003). The discussion above suggests that the structures A, B, and the E2 EELR can 
be identified with a tidal tail driving material to the central region most probably due to a 
past merger event in 3C 120. Most of the analyzed continuum dominated structures shows 
a rather young stellar population, indicating a recent (but not ongoing) star formation rate. 
This result agrees with the picture of a post-merging event. Similar results have been found 
in the study of the optical colors of host galaxies AGNs (Jahnke, Kuhlbrodt, & Wisotzki 
2004; Sanchez et al. 2004b) . 

The poor spatial coverage of the third identified EELR, E3, prevents us of drawing a 
clear picture of its origin. In fact, this region corresponds (sec section §5.1) to a tail extension 
of a larger structure already reported in this object (Soubeyran et al. 1989). Although we 
cannot be conclusive with our current data, the kinematical results most probably indicate 
that E3 is associated with an infiow/outfiow of gas to/from the circumnuclear region. 

It has long been suggested that strong interactions and galaxy mergings may (re-)ignite 
nuclear activity (e.g. Sanders et al. 1988). Galaxy interactions can produce the loss of 
momentum required to allow the infall of gas towards the nuclear regions, gas that would 
feed the AGN. Many authors found that AGN hosts show distorted morphologies, reminis- 
cent of past merging events (e.g. McLeod & Rieke 1994a,b; Bahcall et al. 1997; Sanchez & 
Gonzalez- Serrano 2003). They are found in environments with high probability of experienc- 
ing interactions (e.g. Sanchez & Gonzalez-Serrano 2002). And their hosts, mostly early-type 
galaxies, present anomalous bhic colors (e.g. Sanchez ct al. 2004a). All together it may 
indicate that, if not all, at least a family of AGNs is generated by the merging/interaction 
between galaxies (Canalizo & Stockton 2001). A merging event alone id not enough to 
generate an AGN. The presence of a massive black-hole in the progenitor galaxy is a basic 
requerement. Massive black-holes are only found in bulge-dominated massive galaxies, due 
to the black- hole/bulge mass relation (e.g. Magorrian et al. 1998; Kormendy & Gebhardt 
2001). Indeed, recent results (Sanchez et al., in prep.) show that the fraction of AGNs 
increases in galaxy mergers between two large galaxies or a large with an small galaxy. 

Our current results agree with this sceneario. 3C 120 is a bulge dominated galaxy which 
has, most probably, experienced a merging event with a less massive galaxy. That galaxy was 
completely disrupted in the merging process, falling in parts which produce, most probably, 
many of the observed structures (Moles et al. 1988). This may explain the different stellar 
populations of the structures and the average stellar population in the object (§4.2) as it 
has been previously reported for other galaxies (Prieur 1990). A subtantial fraction of its 
gas has been channelled towards the inner regions, following the detected arm-like structure, 
and concentrating in the E2 region. 
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8. Summciry 

We obtained integral field optical spectroscopy of the Seyfert 1 radio galaxy 3C 120. 
The homogeneous data, excellent spatial and spectral coverage, and good spatial and spectral 
resolution make this atlas a useful tool for studying 3C 120 in the optical. These IFS optical 
data were combined with high resolution HST imaging. The analysis of these data suggests 
that a Seyfert 1 nucleus at the center, an early type galaxy, and several structures formed as 
a consequence of a merging process in the past history constitute the radio galaxy 3 120. At 
least one of these structures is identified with an infiow, which is feeding the central engine. 
A radio-jet is escaping from the center, perturbing the gas on its path. 

The main results from the analysis of this dataset were the following: 

1. Several continuum-dominated structures were detected in 3C 120, which do not follow 
the mean distribution (bulge dominated) of the stellar component. Some of these 
structures were shells in the central Kpc of the object, which may indicate a past 
merging process. The colors of the different components were not compatible with a 
single stellar population in this galaxy. 

2. Three emission line structures were identified (Ei, E2 and E3) which were not associated 
with the general behavior of the galaxy. These gaseous structures presented a high level 
of ionization. The origin of structure Ei was the interaction between the intergalactic 
medium and the radio-jet emerging from the nucleus. 

2. The spectra of the nucleus and the host galaxy were decoupled, obtaining a data cube 
of residuals. Prom that residual data cube, we extracted and analyzed the spectra of 
the different structures in 3C 120. 

3. The velocity field indicated a rotational component plus several kinematical perturba- 
tions associated with the identified emission structures. 

4. The continuum-dominated structures A, B, the Ss shell and the EELR E2 seem to be 
physically associated, belonging to the inner-most part of a tidal-tail, remanent of a 
past merging event. This tail has channeled gas into the inner regions, in an inflow, 
that has generated E2. 

The 4.2-m WiUiam Herschel Telescope is operated by the Isaac Newton Group at the 
Observatorio de Roque de los Muchachos of the Instituto de Astrofisica de Canarias. We 
thank all the staff at the Observatory for their kind support. This project is part of the 
Euro3D RTN on IFS, funded by the EC under contract No. HPRN-CT-2002-00305. This 
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A. Atlas of spectra 

In this appendix we present the spatial distribution of the brightest emission hne profiles 
as spectra diagrams. These diagrams represent line profiles in a short wavelength range of 
the individual spectra at each point (fiber). The spectra at each location are autoscaled to 
show the profile shape (lines nearer to the optical nucleus are brighter than those farther 
out). Figure 13 and 14 show the spectra diagrams corresponding to H/3+[OIII]AA4959, 5007, 
and Ha+[NII]AA6548,6584+[SII]AA6716,6730, respectively. A considerable blending of the 
emission lines with the wings of the broad component of the Balmer lines can be clearly 
appreciated in these diagrams. 

B. Decoupling the spectra of different components using IFS 

A new technique to decouple the spectra of different components in an object has been 
developed, based on well known techniques applied in the decoupling of different components 
in images. Two different methods have been developed. A first one is based on a 2D modeling 
of the objects using object templates, while the second one is based on an isophotal analysis 
of 2D images. We explain here how both methods have been extended to the analysis of IFS 
datacubes, and the advantages of each one. 

B.l. 3D modeling 

IFS combines the characteristics of imaging and spectroscopy. In general terms, IFS 
data consist of a series of spectra obtained at a discrete number of positions in the sky. 
In the case of INTEGRAL data in the SB2 configuration, they consist of the 209 spectra 
obtained in the same number of positions in the sky. Newer instruments, like PMAS (Roth 
et al. 2000) or VIMOS (LeFevre et al. 2003), sample the sky in a continuous way, by coupling 
the fibres to lens-arrays. In these latter cases a IFS datacube can be understood as a stack of 
narrow-band images. More generally, at each wavelength we have a 2D discrete distribution 
of fluxes, that can be transformed to a regular grid distribution by an interpolation routine. 

It is a quite common situation that an observed object is the combination of different 
components: e.g., an QSO lens, a crowded-field of stars and an AGN-|-host system. Different 
techniques have been developed to decouple the different components of an image, and 
to derive their photometry. A wide extended technique is to fit the images with a 2D 
model, including a template for each different component. This technique is implemented 
in GALFIT (Peng et al. 2002), a program for modeling multi-components in images. This 
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program has been used in the decouphng of nuclear and host components (Sanchez et al. 
2004b) and the decouphng of different components in a QSO lens (Wisotzki et al. 2004). 
GALFIT fits the images with a multi-component model, producing a template image and a 
residual image. This residual image can be used to detect substructures in 2D images. 

A natural extension of the 2D modehng technique to IFS is to split the datacube in a 
set of narrow-band images of the width of the spectral pixel (the so-called monochromatic 
images), and treat them as individual images. This technique has been used successfully for 
the deblending of QSO lenses with IFS (Wisotzki et al. 2003). In the particular case of our 
INTEGRAL data, wc first transformed the discrete set of spectra in a regular grid datacube 
using an interpolation routine implemented in E3D (Sanchez 2004). The final datacube has 
a O^'SxO'/S pixel size, i.e., l/3th of the original diameter of the fibers. 

The 2D image modeling of the nucleus and the host for each monochromatic image was 
performed using GALFIT. The 2D model comprises a narrow Gaussian function (to model 
the nucleus) and a de Vaucouleurs law (to model the galaxy), both convolved with a PSF. The 
PSF at each wavelength was obtained from a cahbration star datacube, observed just before 
the object. In a first attempt all the parameters of the model were fitted: the intensities of 
the nucleus and the host, the centroid of the object, the effective radius, the position angle 
and the ellipticity of the host galaxy. As already noticed by Wisotzki et al. (2003), this 
first modeling does not produce acceptable results. The large number of parameters, the 
degeneracy between some of them (e.g., the total intensity and the scale- length), and the 
limited sampling of our data impose limitations to the quality of the recovered spectra. 

However, Wisotzki et al. (2003) demonstrated that it is possible to increase the quality 
of the modeling imposing certain conditions in the structural parameters. The centroid of 
the object varies through the wavelengths due to differential atmospheric refraction (DAR) 
(Filippenko 1982). Therefore, whenever there is a clear peak in the object (like in our 
case) the centroid should vary smoothly with the wavelength. Figure 15 shows the centroid 
coordinates as a function of the wavelength resulting from the 3D modeling. Despite of the 
fiuctuations, there is a clear trend. The sohd-lines show the result of a fitting a polynomial 
function of order 5. In a second run of the 3D modehng we imposed the centroid to be fixed to 
the result from this polynomial fitting. Similar conditions can be imposed for the ellipticity, 
effective radius and position angle of the host galaxy. By forcing the structural parameters to 
vary smoothly with the wavelength and fixing them in the 3D modeling process the number 
of parameters is reduced to the intensity of the different components. Thus the accuracy of 
the modeling is increased (Wisotzki et al. 2003, 2004). 

The number of parameters can be reduced too if wc derive the structural parameters 
from an external source of information. As we explained above we have good quality HST 
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images of the object in different bands, that have been already analyzed using GALFIT 
(§4.1). We fixed the structural parameters (effective radius, ellipticity and position angle) in 
our second 3D modeling of the values derived by the 2D modeling of the HST images (Table 
2). The centroids were fixed to the values derived using the polynomial function shown in 
Fig. 15. Figure 16 illustrates the process. The top-left panel shows the decoupled spectra of 
the nucleus and the host derived from the 3D modeling when all the parameters are fitted 
freely. The spectra derived from the first run are noisy and can only be considered as a 
very rough approximation. The top-right panel shows the spectra once fixed the centroids 
to the result of the polynomial function fitting shown in Fig. 15. There is a clear increase 
of the quahty of the derived spectra. The bottom-left panel shows the spectra when, in 
addition to the centroids, the structural parameters of the host galaxy (ellipticity, position 
angle and cfi^cctivc radius) arc fix to the result from a polynomial function fitting similar to 
the one used for the centroids. This is the best possible result if we only had the information 
from the IFS. The bottom-right panel shows the final result, obtained when the structural 
parameters are fixed to the values derived from the 2D modeling of the HST images. The 
derived spectra are remarkable similar to those of the previous case. However, in the previous 
case it is possible to identify some fake absorption features in the spectrum of the nucleus 
coincident with the narrow emission fines (e.g., [OIII]A5007). We cross-checked the accuracy 
of the final decoupled spectra by comparing the fiuxes derived with that obtained by the 2D 
modefing of the HST images. Both fiuxes agree within ~20%. 

As we explained above GALFIT provides us with a residual image at each wavelength 
once subtracted the 2D model from each monochromatic image. We stored and stacked 
those residual images to create a residual datacube. The residual datacube is used (1) to 
cross-check the accuracy of the fitting, i.e., no evident fake residuals are found, and (2) to 
study the substructures in the object. Figure 17, left panel, shows a narrow-band image 
centred in the continuum adjacent to the [0111] emission line (5204-5246A) from the residual 
datacube (contour-plots), together with the residual image obtained by a 2D modcUng of 
the HST F555W-band image (grayscale). A simple comparison of both images shows that 
the residuals are remarkable similar outside ~2" from the nucleus. However, there is a ring 
structure at ~1'.'5, ~100 fainter than the nucleus, expecially bright in the S5 of the nucleus. 
This structure coincides somehow with the south-east shell. However it is worrisome that it 
is fainter in the north-west, where the Sn shell is brighter. Taking into account the original 
size of the fibers (~0'.'9) and our limited determination of the PSF, we think that this ring 
substructure is an artifact of the 2D modeling. Although it limits the use of the residual 
datacube to study the properties of the substructures, it does not significantly affect the 
accuracy of the decoupled spectra. 



-27- 



B.2. 3D surface brightness analysis 

We developed a second method to create a 3D template of the object that would not 
depend on the PSF determination and without assuming any certain law to describe the 
intensity profile. This method is an extension of the isophotal surface brightness analysis 
technique (Jedrzejewski 1987). It looks for the ellipse parameters that better represent the 
isophotal shape by a Fourier analysis of the surface brightness at a certain radius as a function 
of the eccentric anomaly. The output is a surface brightness profile of the object (i.e., the 
radial distribution of the intensity), and the radial distributions of the centroid, position 
angle and ellipticity (understanding radius as semi- major axis of the ellipse). A template of 
the object is created using these output parameters. Then, a residual of the image is obtained 
by subtracting this template from the original image. This technique has been extensively 
used for the detection of substructures in galaxies for decades (e.g. Perez-Fournon et al. 1986; 
Hjorth et al. 1995). 

We extended this method to analyze substructures in IFS data. For each monochromatic 
image we performed an isophotal surface brightness analysis using our own coded routines. 
The average centroids, ellipticities, position angles, and the integrated fiuxes obtained from 
this analysis are stored for each wavelength. The residual images obtained are combined in a 
residual datacube. Like in the previous method the quality of the model template can be in- 
creased assuming that the structural parameters change smoothly with the wavelength. The 
structural parameters for each wavelength were fixed to the values derived by a polynomial 
function fitting. In the case of an isophotal analysis, once fixed the structural parameters, 
the derived model for each monochromatic image is just estimated obtaining the average 
fiux at a certain radius along the eccentric anomaly. I.e., it is not a fit anymore, but just a 
direct measurement. 

Figure 18 illustrates the results of this process. It shows, in both panels, the integrated 
spectrum of 3C 120 over the field-of-view of INTEGRAL (black line), together with the 
recovered spectrum from the 3D template (orange line), and the integrated spectrum of the 
residuals. The left panel shows the result from the first iteration where all the structural 
parameters are derived by the described method. In general terms the spectrum of the model 
matches well the original spectrum, despite some wavelengths where the program is unable 
to derive the correct parameters. The right panel shows the result from the second iteration 
with all the structural parameters fixed, as described above. In this case the agreement 
between the recovered and the original spectrum is remarkable good, with no flux deviation 
higher than a few percent at any wavelength. 

As expected this method produces a better quality residual datacube than the 3D mod- 
eling. The basic reason is that we reduced the number of assumptions about the shape of 
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the object and we transformed a fit to a direct measurement. Figure 17, rigfit panel, sfiows 
a similar narrow-band image of the continuum adjacent to [OIII] from the residual datacube 
than the one presented in the left panel, but obtained with the 3D surface brightness analy- 
sis. The grayscale shows once more the residual from the 2D modeling of the F555W-band 
image. The residuals from the narrow-band image and the HST data match better than that 
of the 2D modeling. No trace of the ring structure is present, which reinforce our idea that 
it is an spurious result due to a limited sampling and inaccuracies in the PSF determination. 

The two developed techniques to decoupling spectra of different components from IFS 
data give successful results. The first method could be useful to dcblcnd spectra, such in the 
case of spectra from different components. The second procedure is more suitable to study 
the residual structures once the spectrum of the main component is subtracted. 
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Table 1: Summary of ]:)r()])erties of the HST images 



Band 


Exp.Time 


3cr lim. mag 


mag 


mag 




(s) 


(mag/arcsec^) 




NED 


F555W (V) 


2040 


25.3 


14.94 


14.7 


F547M (V) 


2300 


25.1 


14.35 




F675W (R) 


2250 


25.3 


13.74 


13.9 


F814W (I) 


2236 


24.9 


(1) 


12.9 


(1) Peak saturated. 



Table 2: Summary of the result from the 2D modeling 

Band n mnuc mhost a/b PA(") 

F555W 4.7 16.1 15.4 1.9 kpc 0.87 -16.8 

F547M 4.5 16.9 14.5 2.9 kpc 0.86 -18.6 

F675W 6.7 16.1 13.9 2.8 kpc 0.73 -19.1 

F814W — — 13.2 0.94 -18.5 
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Table 3: Properties of the emission lines detected in the nucleus and host galaxy of 3C 120 



Name 



Nucleus 



HOST 



Line 




Ai 




o-A 


Flux2 








4250, 


.33 


24, 


.12 


61, 


.25 


260, 


.78 


H5BL1 


4249, 


.15 


24, 


,12 


319, 


,60 


583, 


,19 


H(5BL2 


4254, 


.50 


24, 


.12 


1754, 


.62 


516, 


.10 


H7 


4498, 


.42 


24, 


.12 


35, 


.08 


268, 


.68 


H7BLI 


4497, 


.24 


24. 


.12 


885, 


.98 


539, 


.97 


H7BL2 


4502, 


.58 


24, 


,12 


790, 


,73 


422, 


,00 


[OIII]A4363 


4521, 


.98 


24, 


.12 


60, 


.99 


119, 


.60 


HeIIA4686BLl 


4846, 


.66 


24, 


,12 


192, 


.56 


100, 


.54 


HeIIA4686BL2 


4876, 


.41 


24, 


,12 


270, 


,39 


141, 


,18 


H/3 


5024, 


.82 


0, 


,19 


75, 


,54 


479, 


,99 


H/3BL1 


5028, 


.23 


1. 


,76 


1473, 


.31 


754, 


.88 


H/3BL2 


5058, 


.85 


53, 


.77 


1676, 


.73 


788, 


.82 


[OIII]A4959 


5125, 


.77 


0, 


,19 


163, 


,17 


35, 


,06 


[OIII]A5007 


5175, 


.36 


0. 


,19 


489, 


.58 


105, 


.20 


HeIA5876BLl 


6075, 


.86 


11. 


.62 


332, 


.05 


112, 


.13 


HeIA5876BL2 


6105, 


.86 


11. 


.62 


465, 


.42 


157, 


.05 


[OI]A6300.3 


6518, 


.52 


40, 


.42 


20, 


,34 


33, 


,63 


[OI]A6363.8 


6584, 


.19 


40, 


.42 


30, 


.34 


34, 


.81 


[NII]A6548 


6784, 


.48 


2. 


,08 


9. 


,46 


57, 


.34 


Ha 


6799, 


.71 


2. 


,08 


232, 


.62 


513, 


.60 


HaBLl 


6782, 


.94 


0. 


,58 


3570, 


.51 


574, 


,07 


HaBL2 


6797, 


.59 


1. 


,15 


5985, 


.15 


446, 


.94 


[NII]A6583 


6821, 


.21 


2. 


,08 


28, 


.40 


172, 


.13 


[OII]A3727 


38G7, 


,81 


10, 


,22 


477, 


,51 


60, 


,96 


[NeIII]A3869 


4010, 


.57 


16, 


,32 


332, 


,54 


48, 


,24 


[NeIII]A3967 


4117, 


.10 


16, 


.32 


179, 


.77 


34, 


,32 


[SII]A4072 


4223, 


.39 


122, 


.12 


110, 


.88 


37, 


.32 


H7 


4494, 


.84 


389, 


,92 


54, 


,84 


28, 


,56 


[OIII]A4364 


4518, 


.40 


389, 


.92 


56, 


.88 


28, 


.32 


HeIIA4686 


4851, 


.17 


143, 


.31 


101, 


.76 


33, 


.48 


H/3 


5030, 


.04 


0. 


,82 


218, 


.05 


39, 


.72 


[OIII]A4959 


5131 


.00 


0, 


,82 


1164, 


.08 


161, 


,37 


[OIII]A5007 


5180, 


.58 


0. 


,82 


3495, 


.62 


484, 


.59 


HeIA5876 


6072, 


.99 


23, 


.27 


45, 


.60 


47, 


.40 


[OI]A6300 


6514, 


.63 


194, 


,46 


122, 


,40 


60, 


,48 


[OI]A6364 


6580, 


.30 


194, 


,46 


58, 


,08 


55, 


,20 


[NII]A6548 


6771, 


.69 


2. 


,31 


259, 


.45 


46, 


.92 


Ha 


6786, 


.91 


2. 


,31 


2097, 


.63 


255, 


.97 


[NII]A6583 


6808, 


.41 


2, 


,31 


778, 


,37 


140, 


,77 


[SII]A6716 


6946, 


.13 


32, 


,68 


368, 


,18 


96, 


,00 


[SII]A6731 


6961, 


.00 


32, 


,68 


248, 


,29 


89, 


,52 


Arm 


7377, 


.62 


291, 


.60 


108, 


.48 


69, 


.00 



(1) wavelength in A, (2) flux in lO'^^ 



erg s cm units, (3) the narrow-emission lines had a velocity dispersion of ~700 km s , while the 



broad-emission lines have to component, BLl and BL2, with a velocity dispersion of - 

^8706 km s~^, respectively. 



.^2188 km s-^and 



-36- 



Table 4: Properties of the emission lines detected in the EELRs of 3C 120 
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Tablc 5: Properties of the gaseoTis emission of the different components in 3C f 20 
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Figure 1: F555W-band image of 3C 120 taken witli tlie WFPC at tlie HST. Tiie spatial 
distribution of the science fibers at the focal plane of the 4.2m WHT has been overlaid. 
North and East are up and left, respectively as usual. The orientation is the same in all 

figures in this paper. 
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Figure 2: Nuclear spectrum in the full spectral range. This spectrum corresponds to a 
hexagonal aperture of 1'.'6 in radius centred on the optical nucleus. 
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Figure 3: (A) Two-dimensional distribution of the continuum emission of 3C 120 derived 
from the IFS data by integrating the signal in the spectral range 5100-6100 A, similar to 
that range covered by the IZ-band filter. (B) F555W-band image of the same central region 
of 3C 120 obtained with the HST/PC (O'.'046/pixel). Different continuum dominated 
structures already detected from ground base observations are labeled using the 
nomenclature from Soubeyran et al. (1989). The contour-plot shows the radio map at 4885 
MHz Walker (1997). (C) F-band residual image, once subtracted the galaxy template 
derived by a 2D modeling from the y-band image shown in panel (A). (D) F555W-band 
residual image derived by subtracting a galaxy template obtained by a 2D modeling from 

the F555W-band image shown in panel (B). 
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Figure 4: (A) The V — I color image derived from the V and J maps recovered from the 

IFS data. (B) The V — I color image obtained from the F555W and F675W-band HST 
images. The image has been smoothed with a 3x3 pixel median kernel. (C) The residual V 
map (close similar to an ionized gas map) derived subtracting a continuum recovered from 

the IFS data (5600-5700 A) avoiding the contribution of any bright emission line and 
scaled. (D) The residual image obtained by the subtraction of the F547M-band image from 
the F555W-band image, once scaled by the differences of the photometric zero-points. The 
contour-plot corresponds to the F555W-band residual image shown in figure 3D. Labels 
indicate the previously detected EELRs (Soubeyran et al. 1989; Sanchez et al. 2004a). 
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Figure 5: V — R colors as a function of the R — I colors of the host galaxy and the 
different structures found in 3C 120. The discontinuous line shows the location of a single 
stellar population model at different ages Bruzual & Chariot (2003). The values indicate 

the logarithm of the age in Gyrs. The dotted line shows the average colors of different 
galaxy types (E, Sab, Sbc, Scd and Irr), obtained from Fukugita et al. (1995). The arrow 
shows the effect of dust in the colors for an absorption of Ay ~1.5 mags. Errors of the 
individual colors are 0.1 and 0.26 mag in V-R and R-I, respectively. 
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Figure 6: Integrated intensity maps of (A) [OIII]A5007 and (B) Ha emission lines derived 
by a Gaussian fitting. In the circumnuclear region, we have included a second broad 
Gaussian to Ha. Dividing line intensity maps, we have derive: (C) the ionization map 
[0III]/H6eta; (D) [Nil] /Ha 2D-distribution; (E) the extinction Ha/H/3 map; and (F) the 
electronic density distribution [SII]A6716/[SII]A6730. 
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Figure 7: (A) Contour-plot of a narrow-band image centred on tlie continuum adjacent to 
the [OIII] emission line (5204-5246A) extracted from the residual data cube, together with 
a greyscale of the residual from the 2D modeling of the F555W-band image (Fig. 3D). (B) 
Contour-plot of a narrow-band image centred on the [OIII]A5007 emission line 
(5170-5200A) extracted from the residual data cube, together with a greyscale of the 
F555W-F547M image shown in Fig. 3D. The different continuum dominated (left) and 
emission-line dominated (right) substructures are labelled with their corresponding names. 
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Figure 8: (A) Spectrum of the nucleus of 3C 120, once decontaminated from the 
contribution of the host galaxy. (B) Spectrum of the host galaxy (black line) together with 

the spectrum of synthetic models for single stellar populations of 16 Gyr (red line), 0.01 
Gyr (blue line), and a mix of both (orange line), adding (solid hne) and not adding (dashed 
line) a power-law continuum to simulate scattered light from the nucleus. 
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Figure 9: (A) Spectra of the different EELRs detected in the residual data cube. (B) 
Spectra of the different continuum-dominated structures detected in the residual data cube 
(points) and a median filtered version of these data (solid lines, solid squares), with a 300A 

width. 



-46- 



-| 1 1 r 



1 n 1 r 



J I I L 



T 1 1 r 




J I L 



J I I L 



d 



ID 



(s'h/^oo9y[iiioJ^§oi 



Figure 10: [OIII]A5007/H/3 vs. [NII]A6583/Ha intensity ratios for the different 
components of 3C 120. Solid curve divides AGNs from Hn region-like objects. 
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Figure 11: Ionized gas velocity fields of 3C 120 derived from the Gaussian fit to: (A) 
[OIII], isovelocity lines extend from 10200 to 10600 in steps of 50 km s~^; (B) Ha where 
isovelocity lines extend from 10000 to 10600 in steps of 50 km s— 1. (C) As (A) with the 
[OIII] intensity of Ei, E2, E3 in black contours. (D) The velocity distribution of ionized gas 
structures (Ei, E2, E3). Isovelocity lines extend from 9600 to 10300 in steps of 100 km 
s^^(E) The simplest rotational model obtained considering an inclination of 40° and 
PA=-40°. (F) The residual velocity map after the subtraction of the rotational model (E) 
from the ionized gas velocity field (A). Isovelocity lines corresponds to kinematics of 

structures Ei, E2, E3 as in (D). 
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Figure 12: The solid line shows the Ha velocity distribution of 3C 120 along a position 
angle of —40° from south-east (left) to north-west (right). The perturbation associated 
with E2 is seen as a peak at approximately —5". The dashed-line shows the distribution 
obtained assuming that the north-west portion of the velocity distribution describes the 
behavior of the rotating component of the galaxy, and modeling the south-east by a 

symmetrical distribution. 
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Figure 13: Two-dimensional spectroscopic diagram of H/?+[OIII]AA4959, 5007 emission 
line profiles. Plotted spectral range: 4910-5300 A. 
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Figure 14: Two-dimensional distribution of emission line profiles of 
Hq;+[NII]AA6548,6584+[SII]AA6716,6730. The spectral range plotted is 6650-7050 A. 
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Figure 15: Centroid coordinates, x and y, in pixels ((/.'S /pixel), as a function of the 
wavelength as derived from the first run of the 3D modeling. The solid-lines show the 
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result of the fitting with a polynomial function. 
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Figure 16: Top-left panel: Decoupled spectra of the nucleus and the host of 3C 120, 
together with the total spectrum of the object, derived from the 3D modeling when all the 

parameters are fitted freely. Top-right panel: Similar spectra derived from the modeling 
when the centroid is fixed to the result of the polynomial function fitting shown in Fig. 15. 
Bottom-left panel: Similar spectra derived when in addition to the centroids, the ellipticity, 
position angle and effective radius of the host galaxy are fix to the result of a polynomial 

function fitting as a function of the wavelength. Bottom-right panel: Final decoupled 
spectra obtained when the structural parameters of the host galaxy are fixed to the values 
derived from the 2D modeling of the HST images. 
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Figure 17: Left panel: Contour-plot of a narrow-band image centred on the continuum 
adjacent to the [OIII] emission line (5204-5246A) from the residual datacube obtained by 

the 3D modeling. Right panel: Contour-plot of a narrow-band image centred on the 
continuum adjacent to the [OIII] emission line (5204-5246A) from the residual datacube 
obtained by the 3D surface brightness analysis. In both panels the grayscale shows the 
residual from the 2D modeling of the F555W-band image (Fig. 3D) 
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Figure 18: Left panel: Integrated spectrum of 3C 120 (black-line), together with the 
spectrum of the model (orange) and the residuals (green), obtained using the 3D surface 
brightness analysis. Right panel: Similar spectra than the ones presented in the left panel 
obtained once the structural parameters are fixed to the result of a polynomial function 

fitting as a function of the wavelength. 



